Vortex phase diagram of Ba(Feo.93Coo.o7)2As2 single crystals 
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Detailed measurements of the global and local electromagnetic properties of Ba(Feo.93Coo.o7)2As2 
single crystals are reported. Analysis of the irreversible magnetic response provides strong evidence 
for similar vortex physics in this Fe-based pnictide superconductor to the high- Tc cuprates, such as 
Y-Ba-Cu-0 or Nd-Ce-Cu-O. In particular, we have found a nonmonotonic "fishtail" magnetization 
in M {H, T = const) loops and its signature is also present in AI {H = const, T) scans. The super- 
current density is evaluated by using several techniques, including direct transport measurements. 
At 5 K we estimate its value to be a moderate j « 2.6 ± 0.2 x 10^ A/cm^. Analysis of the magnetic 
relaxation is consistent with the collective pinning and creep models (weak pinning and fast creep) 
and suggests a crossover from the collective to the plastic creep regime in fields exceeding the value 
corresponding to the maximum in fishtail magnetization. 

PACS numbers: 74.25.Sv, 74.25.Ha, 74.25.Qt, 74.25.Dw, 74.25.-q 



I. INTRODUCTION 

The behavior of type-II superconductors in applied 
magnetic fields has been a subject of ongoing inter- 
est due to obvious applications as well as deep connec- 
tions to the microscopic properties of superconductors 
and mechanisms of superconductivity (for some reviews, 
see Refs. [lIM's'.'d'.'s'.'o'.T]). The situation is complicated 
in type-II superconductors because magnetic fields pen- 
etrate the sample in the form of Abrikosov vortices^iS 
whose collective behavior determines the macroscopic 
electromagnetic response. The interaction between vor- 
tices and local modulations of the superconducting prop- 
erties leads to pinning, which together with thermal fluc- 
tuations results in spatio-temporal variations of the vor- 
tices (magnetic induction) and therefore the irreversible 
magnetic properties. Strong thermal fluctuations and 
weak pinning lead to complex H — T diagrams with vari- 
ous crossovers in the static and dynamic electromagnetic 
response of vortex matter. 

In the case of layered superconductors, such as the 
high-Tc cuprates and ET organics, a non-monotonic de- 
pendence in the field dependent magnetization is ob- 
served (see Refs. Pllillijlil). This is called the "fish- 
tail" or peak-effect. The latter term is ambiguous since 
a different nonmonotonic behavior, usually closer to the 
second critical field, is observed in low- Tc superconduct- 
ing alloysiii^. Examples of this can be found in materi- 
als such as NbSe2^ and neutron-irradiated VaSii^ and 
MgB2.^^ Another prominent feature of high- Tc supercon- 
ductors is a very large rate of magnetic relaxation (giant 
flux creep, reported by Yeshurun and Malozemoff^^) and 
its complicated time, magnetic field, and temperature 
dependence. ^'^'^ While the mechanisms of the magnetic 
relaxation and fishtail feature in the cuprates are still 
being discussed, several self-consistent microscopic nu- 
merical simulations and phenomenological models have 
been proposed. Some of these have been able to explain 
the observed effects in many classes of unconventional 



superconductors. Magnetic relaxation with an empha- 
sis on high- Tc superconductors is reviewed by Yeshurun, 
Malozemoff, and Shaulov.^ A detailed analysis of the pin- 
ning strength and the barriers for magnetic relaxation in 
the case of weak collective pinning and creep is given by 
Blatter et al^ and also by Brandt.^ 

In disordered granular and polycrystalline materials 
much of the response comes from extrinsic factors, such 
as granularity and significant disruption of supercurrent 
flow between the grains. To study intrinsic magnetic 
properties that can be related to the mechanism of super- 
conductivity, such as anisotropy and basic length scales 
(coherence length, ^, and London penetration depth. A), 
one must turn to single crystals. A vast amount of litera- 
ture on this subject clearly shows that different classes of 
superconductors exhibit unique static and dynamic mag- 
netic responses. This allows us to discuss the similarities 
and differences in vortex behavior and to look for the 
connections between different classes of materials. Vor- 
tex phase diagrams that show various transitions, such 
as changes in vortex lattice symmetry, dimensionality or 
mobility. Pinning strength and magnetic dynamics are 
particularly useful for such comparisons and determin- 
ing applications. Recent examples are CeCoIng^ and 
MgB2^. 

So far, only a limited amount of information regard- 
ing the vortex behavior in single crystals of Fe-based 
pnictide superconductors is available. In polycrystalline 
samples, the magnetic behavior is significantly affected 
by extrinsic factors, such as grain morphology, surface 
roughness, and inter-grain voids and interfaces. In these 
materials, it has been found that some of the grains are 
not even superconducting^^. Although the superconduct- 
ing granularity can be estimated from magneto optical 
imaging, and some important parameters can still be ex- 
tracted^iS?^ the polycrystalline structure docs compli- 
cate the analysis. For example, a signature of nonmono- 
tonic M{H) has been reported in polycrystalline wires of 
SmFeAsOo.8Fo.2^'* and could be related to the fishtail fea- 



ture of this work, but it could also be due to the physics 
of inhomogeneous multi-phase materials. Furthermore, 
severe effects of inter-grain weak-links did not allow for 
quantitative analysis of the intrinsic properties. 

The first measurements of magnetization and relax- 
ation in individual single crystals of the oxypnictide su- 
perconductor NdFeAsOi^ajFa; have revealed the vortex 
physics to be surprisingly reminiscent of cuprate super- 
conductors, in particular the existence of a large mag- 
netic relaxation rate and relatively weak pinning. Later 
torque measurements on the similar Sm-based system 
have lead to the conclusion that the electromagnetic 
anisotropy is large and quite temperature dependent. 
More recent high-field torque measurements on similar 
crystals have shown a much smaller anisotropy^- The su- 
perconducting crystals based on the parent RFeAsO com- 
pound (1111 system, where R is a rare earth element), are 
small and difficult to obtain (via high-pressure synthesis) 
and study. On the other hand, large crystals based on the 
related oxygen-free parent compound AFe2As2 (122 sys- 
tem, where A is an alkaline element), can be grown in flux 
at ambient pressure (see, e.g., A—Ba^, Sr— and Ca3i). 
Both A and Fe sites can be doped to achieve supercon- 
ductivity with holes or electrons as carriers, respectively. 
For example, among other results, detailed thermody- 
namic measurements of the Bai_2;K2:Fe2As2 compound 
have shown a very low anisotropy of the second critical 
field, ~ 321i2£. The work performed on superconducting 
Ba(Fei_a;Coa;)As2 has shown that Co doping does not 
introduce significant additional scattering'^^. The fishtail 
magnetization and weak anisotropy has also been stud- 
ied. High-field properties''^ as well as penetration depth 
studies'^* of single crystals of this compound will be re- 
ported elsewhere. It is this system that we chose for the 
present study for its superior superconducting properties. 

In this paper we use global and local magnetic prop- 
erties as well as direct electro-transport measurements 
to study the details of vortex pinning and magnetic re- 
laxation, evaluate supercurrent density, and finally con- 
struct the vortex phase diagram in single crystals of 
the recently discovered Fe-based pnictide superconduc- 
tor, Ba(Feo.93Coo.o7)2As2. 

We distinguish between the true critical current den- 
sity, jc, which marks the crossover between flux flow 
and flux creep regimes, and the general supercurrent 
density, j < jc, which is observed in the experiment. 
By definition of the critical current, at j = jc the 
current-dependent barrier for vortex escape from the pin- 
ning potential U{j/jc) = 0. The supercurrent is deter- 
mined by the time-window of the experiment, At, and 
can be estimated from the logarithmic solution^^ of the 
flux creep equation, U{j/jc) = Tin (1 -I- A</<o), where 
to^ 10-5- 10-'^ s, is the macroscopic characteristic time 
that depends on the sample size and shape as well as on 
the microscopic attempt time, tq ~ 10"^°- 10"^' s.-i^i^i 
(We use units where the energy barrier is measured in 
kelvin, so fcs = 1.) It is actually the supercurrent that 
determines the usefulness of a particular superconductor 
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FIG. 1: (Color online) Schematics used in calculations and 
definitions of the sample dimensions and flows of the super- 
currents, (a) the magnetic field is applied perpendicular to 
the ab— plane (parallel to the crystallographic c— axis and 
edge), (b) the magnetic field is applied along the ab— plane 
(shown along the a— edge). 



for current-carrying applications and therefore the mag- 
netic irreversibility should be analyzed both in terms of 
magnetic hysteresis and the rate of magnetic relaxation. 



II. EXPERIMENTAL 

Single crystals of Ba(Feo.93Coo.o7)2As2 were grown out 
of FeAs flux using high temperature solution growth tech- 
niques. More speciflcally, the powders of FeAs and CoAs 
were mixed with Ba in the ratio Ba:FeAs:CoAs=l:3.6:0.4. 
The mixture was placed into an alumina crucible and a 
second catch crucible containing quartz wool was placed 
on top of the growth crucible. Both were sealed in a 
quartz tube under argon and slowly heated to 1180 °C, 
held for 2 hours, and slowly cooled to 1000 °C over a 
period of 36 hours. Once the furnace has reached 1000 
°C the FeAs is decanted from the single crystals using a 
centrifuge. The size of the resulting crystals can be as 
large as 12 x 8 x 1 mm'^. Elemental analysis was per- 
formed using wavelength dispersive X-ray spectroscopy 
in the electron probe microanalyzer of a JEOL JXA- 
8200 Superprobe. It showed the actual Co concentration 
is Co/(Co-}-Fe)=7.4%, so the more accurate formula is 
Ba(Feo.926Coo.o74)2As2, but we round it to the second 
digit in the text. During the course of this study several 
samples from the same batch were measured. As shown 
below, they all have a transition temperature off Tc = 22 
K, as determined from dc magnetization (inset of Fig. [2]) 
and from transport measurements. Fig. [71 

For the measurements of total magnetic moment 
(global magnetic measurement), a cuboid-shaped sam- 
ple (see Fig. [3]) of dimensions, 0.28 x 0.7 x 1.26 mm'^, 
was fixed in a gelatine capsule with a small amount of 
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FIG. 2; (Color online) Magnetization loops measured in a 
0.28 X 0.7 X 1.26 mm^ single crystal of Ba(Feo.93Coo.o7)2As2 
along the crystallographic c— axis (along the d— edge) at 
T = 5, 10, 15 and 18 K. Left inset zooms at the higher tem- 
perature region showing the data for T = 15, 16, 17, 18 and 
19 K where the evolution of fishtail magnetization is clearly 
seen. Right inset shows a superconducting transition mea- 
sured after zero-field cooling (zfc) at H = 10 Oe. 



FIG. 3; (Color online) Magnetization measured along the fl- 
edge in, perpendicular to the crystallographic c— axis in a 
Ba(Feo.g3Coo.o7)2As2 single crystal. Left inset: comparison 
between M{H) loops at T = 15 K in two orientations. Right 
inset: a photograph of the measured cuboidal sample with 
dimensions, 0.28 x 0.7 x 1.26 mm^. 



Apiezon grease. The capsule was placed inside of a clear 
plastic straw. To check for possible errors due to me- 
chanical misalignment of the sample assembly, several 
key measurements (such as M [H) loops) were repeated 
three times each time after the sample was removed and 
re-assembled. No noticeable variation in the results was 
found. The magnetization measurements were conducted 
in a Quantum Design MPMS magnetometer. The second 
critical field, Hc2 {T), estimated from the onset of su- 
perconductivity, was measured by using a tunnel-diode 
resonator technique42. 

Samples for transport measurements were cut with a 
razor blade into long thin bars having typical dimensions 
of (2 — 3) X 0.1 x 0.02 mm'^. Contacts were made to the 
samples by soldering silver wires with a silver-based al- 
loy and had negligible resistance (less than 0.1% of the 
sample resistance). This allowed for both two and four 
probe measurements to be taken. To sustain a high cur- 
rent density without thermal shock, samples were glued 
with GE-varnish onto an insulating heat sink substrate 
(LakeShore Cryotronics Inc.). Current wires were ther- 
mally anchored to two silver foil heat sinks (see inset in 
Fig. [5]). Current-voltage characteristics, V(I), were mea- 
sured in a Quantum Design PPMS in constant current 
mode. The highest attainable current density was re- 
stricted to keep Joule heating less than 1 K at Tc. 

Magneto-optical (MO) imaging was performed in a ''He 
optical flow-type cryostat utilizing the Faraday rotation 
of polarized light in a Bi - doped iron-garnet indicator 
film with in- plane magnetization''*^. The spatial resolu- 



tion of the technique is about 3 /im with a sensitivity to 
magnetic fields of about 1 G. The temporal resolution, as 
determined by our image acquisition hardware, is about 
30 msec. In all images, the intensity is proportional to 
the local value of the magnetic induction perpendicular 
to the sample surface. 

III. RESULTS AND DISCUSSION 
A. Supercurrent density 

1. Magnetization 

Let us first define the geometry of the experiment. 
We consider a cuboid-shaped crystal with dimensions 
2d < 2a < 2b, as shown schematically in Fig[l] and one 
of the actual samples, 0.28 x 0.7 x 1.26 mm'^, is shown in 
the inset of Fig [3] In Ba(Feo.93Coo.o7)2As2, the crystal- 
lographic ab— plane has the largest area and is parallel 
to the geometric ab— plane of a cuboid. The smallest 
dimension 2d is the sample thickness. When a magnetic 
field is oriented along the crystallographic c— axis (along 
the d— edge), the measured magnetic moment per unit 
volume (volume magnetization) is denoted as Md. Simi- 
larly, magnetic moments measured along the b— edge and 
along the a— edge are Mi, and Ma, respectively. Due to 
tetragonal symmetry, the induced supercurrents can be 
considered isotropic in the ab— plane (we use the desig- 
nation j for the in-plane supercurrent density) , but may 
be different for the magnetic field oriented perpendicular 
to the c— axis. For the vortex motion crossing the Fe-As 
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FIG. 4: (Color online) Magneto-optical images of a single 
crystal Ba(Feo.93Coo.o7)2As2 (0.053 x 0.24 x 0.55 mms) in 
the remanent state. (Magnetic field was applied along the 
c— axis (perpendicular to the page) above Tc and removed 
after cooling to 5 K). Frames (a), (b) and (c) show images 
for T = 10, 13 and 15 K, respectively. Frame (d) shows a 
three-dimensional Bean oblique wedge where the z— axis is 
intensity corresponding to the frame (a). 



750 
500 



-250 



1 1 1 1 1 


' 1 ' 


B II c-axis 




AB(8.5 K) I f 


\..__- — S.b K 




1 


1 . 1 . 1 . 1 . 



-300 -200 -100 100 200 300 



x(nm) 

FIG. 5: (Color online) Profiles of the magnetic induction mea- 
sured in a single crystal Ba(Feo.93Coo.o7)2As2 at different tem- 
peratures along the line shown in Fig. |4jb). Definition of the 
maximum variation of the magnetic induction, AS is shown 
for profile obtained at T = 8.5 K. 



2. Magneto-optical imaging 



planes, the shielding current density is jj_ and for vor- 
tices moving parallel to these planes we use j\\ . Details 
of the anisotropic response will be published elsewhere4^ 

Figure [2] shows M [H) loops measured at several tem- 
peratures. A noticeable nonmonotonic "fishtail" magne- 
tization develops at elevated temperatures. (It presum- 
ably exists also at the lower temperatures, but is shifted 
to higher magnetic fields beyond the capabilities of our 
setup). The left inset in Figl^] zooms into this temper- 
ature interval to show the details of the fishtail evolu- 
tion with temperature. This behavior is quite similar to 
YBCO single crystals.^liiX 

The situation is quite different when magnetization is 
measured with a magnetic field applied along the crys- 
tallographic ab— plane. Figure [3] shows measurements 
performed with H along the a— edge of the crystal. No 
noticeable fishtail behavior is observed. The inset com- 
pares M{H) loops at T = 15 K for the two orientations 
where the difference is clearly seen. When a magnetic 
field was applied along the longer b— edge, the result 
was the same as in the measurement along the a— edge 
but with a different magnitude due to the difference in 
the cross-section perpendicular to the field. Similar be- 
havior of M {H) in two orientations has been observed 
in several crystals with different aspect ratios varying by 
orders of magnitude, which excludes possible geometric 
effects. 



In order to evaluate the supercurrent density we need 
to examine the structure of the critical state by observing 
the distribution of the magnetic induction in the sample. 
The magneto-optical images of a sample from the same 
batch used in M{H) measurements are shown in FigH) 
The sample dimensions were 0.053 x 0.24 x 0.55 rnvn.^. 
Frames (a), (b), and (c) show the remanent state ob- 
tained at T = 10, 13, and 15 K, respectively. Figure [4] 
(d) shows a 3D plot where the z— axis is the magnetic 
induction. This oblique wedge shape is what is expected 
from the Bean critical state model'*'^''*^. We can therefore 
use this well-known approach to calculate the supercur- 
rent density from the measured magnetization and the 
profiles of magnetic induction across the sample. 

We can also examine the profiles of magnetic induction, 
B{r), across the sample obtained along the line shown in 
Fig.dfb). Figure [H] shows such profiles obtained from the 
measurements at different temperatures. The shape of 
the profile is typical for a thin slab-like sample. There is 
a clear neutral line due to self-fields that are generated 
by the Bean supercurrents. A graphical definition of the 
maximum variation of (r), used later for the evalua- 
tion of the supercurrent density, is shown for the T = 8.5 
K profile, as an example. 



3. Evaluation of the supercurrent density from various 
measurements 

We now extract the supercurrent density from global 
and local magnetic measurements and compare it to di- 
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rect transport measurements that were performed on the 
crystals selected from the same batch. Current-voltage 
characteristics are shown in Fig. [6l In a single crystal 
that was still quite thick for transport measurements, 
we could only induce the transition to the normal state 
in a very limited temperature interval. This was suffi- 
cient to overlap with the magnetization data to check 
whether our indirect calculations of j result in a compa- 
rable magnitude of the supercurrent. Sample dimensions 
were determined in a calibrated optical microscope and 
we estimate the (systematic) errors on the order of 10% 
for the evaluated supercurrent density. Random errors 
were negligible. 

Figure [5] shows current - voltage, V{j), characteristics 
measured at different temperatures. The Inset of Fig. [5] 
shows a photograph of the sample with the contacts. The 
temperature dependence of the resistance measured at 
different excitation currents is shown in Fig. [71 The resis- 
tance, normalized by the room-temperature value, mea- 
sured by a two-probe technique between current leads 
and in the conventional four-probe configuration in the 
full temperature range, is shown in the inset of Fig. [T] 
The data scales perfectly with the sample geometry, as 
the distance between contacts for the voltage readings 
are different. 

In the normal state the V{j) curves remain linear up 
to j Ri 10 kA/cm^. For higher j, self heating effects raise 
the slope of the curve. Below Tc, a notable curvature in 
V{j) develops at zero bias and eventually a broad region 
of zero voltage appears at the lower temperatures indi- 
cating a true superconducting state. With the increase 
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FIG. 6: (Color online) Voltage as a function of current density 
measured at different temperatures in a bar-shaped, 3 x 0.1 x 
0.02 mm'^, single crystal of Ba(Feo.93Coo.o7)2As2 with current 
applied along the ab— plane. Upper inset: A photograph of 
the sample with contacts. Lower inset: The derivative dV/dj 
used to identify the critical current, j. 
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FIG. 7: (Color online) Resistance as a function of temperature 
for different values of a dc transport current in a 3 x 0.1 x 
0.02 mm"^ single crystal of Ba(Feo.g3Coo.o7)2As2. The onset 
of non-zero resistance was used to determine the supercurrent 
density. Inset: normalized resistance in a full-temperature 
range measured by two- and four- probe methods. 



of the current density, the voltage crosses over to a linear 
dependence as it should be in the normal state. From 
these measurements we have determined the value of the 
critical current density as the point at which the deriva- 
tive, dV/dj , is maximal as shown in Fig. [6l Though this 
definition does not represent the true zero resistivity state 
in the vicinity of Tc, from below the transition this cor- 
responds closely to a point of sharp increase in the mea- 
sured voltage and is easy to determine. In order to avoid 
dynamic effect related to sweeping the current, we also 
estimated the current density from R(T) curves measured 
at different values of the bias current. The R{T) mea- 
surements are shown in Fig. [7] and the extracted current 
densities are compared to those obtained from V{I) and 
magnetic measurements in the inset Fig. [H] The overall 
temperature dependence of the supercurrent density ob- 
tained at different temperatures by different techniques 
are summarized in Fig. [51 

To determine the supercurrent from magnetic measure- 
ments, we use the Bean model in which there is a field- 
independent jii^iii Assuming that j is isotropic in the 
crystallographic ab— plane, there are three different cur- 
rent densities depending on the orientation of an external 
magnetic field with respect to the crystal faces. When 
the magnetic field is oriented along the d— edge (crystal- 
lographic c— axis), Abrikosov vortices move in the ab— 
plane and gradients in their density induce the supercur- 
rent, j, given by 

cMd / a\-i 
The other two measured component of volume magneti- 
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FIG. 8: (Color online) Supercurrent density in same-batch 
single crystals Ba(Feo.93Coo.o7)2As2. "Magnetic" current 
density estimated from Eq. [T] at zero applied field by an- 
nealing of the remanent magnetization (filled circles), from 
M{H) loops measured at different temperatures (squares), 
from magneto-optical profiles using Eq. [2] (stars), and from 
direct transport measurements (open circles). The lines show 
fits to j{T) = j(0) (1 - (T/Tcf)" of the dataset shown in 
full circles with indicated exponents and adjusted values. 
Inset zooms at the region of Tc for easier comparison with 
the transport data obtained from Fig. [5] (circles) and Fig. [7] 
(triangles) . 



zation, Ma and Mf,, can be used to calculate j|| and j±. 
A detailed study of the anisotropic properties, including 
supercurrents, will be published elsewhere.'*'' 

Alternatively, the supercurrent can be calculated di- 
rectly from the magnetic induction profiles, (r), mea- 
sured along the line shown in Fig. 2] (b). It is worth 
noting that the total magnetization and magneto-optics 
are two very different measurements performed under dif- 
ferent experimental conditions. The only common part 
of the analysis is the Bean model used for data analysis. 
If this assumption is invalid, the resulting supercurrent 
densities will be very different. The variation of magnetic 
induction, AB — \Bz (0) — B^ (edge)| in the full critical 
state is given hy^ 



= 77 In — + 2 arctan (277) — arctan (77) 

(2) 

where r] = d/a. In the sample used for MO study, 
■q « 0.22 and Eq. [H yields j « 194AB. Note that if 
instead of using Eq. [5] one would use the straightfor- 
ward Bean model for a semi-infinite slab (we still mea- 
sure the z— component of B on the surface), then the 
estimated current density would be significantly lower, 
j = c/(27ra)AB « 133AB. 



We first analyze the supercurrent density for zero ap- 
plied field. Figure [5] shows the supercurrent density esti- 
mated from magnetization, magneto-optical, and direct 
transport measurements. From magnetization, the cur- 
rent was obtained from two types of measurements. First, 
a full remanent state was induced at 5 K by cooling in a 5 
T magnetic field and then turning the field off. Then, this 
remanent state was slowly warmed up, allowing the tem- 
perature to stabilize before the measurement was taken. 
This convenient method is relatively quick and produces 
a curve with many data points. However, it is unclear 
whether we probe the original critical state relaxed only 
during the time window At, characteristic of our magne- 
tometer, or the supercurrent has relaxed more and we are 
probing a deeply relaxed state. Therefore, we have used 
M{H = 0) values directly from the magnetization loops 
measured at different temperatures and plot the result- 
ing i as open squares in FiglH The data simply falls on 
top of the curve obtained from the annealing of the re- 
manent state (filled circles). The explanation is that the 
change of temperature occurs during a relatively short 
period of time, comparable to the time of setting and sta- 
bilizing the magnetic field in the superconducting mag- 
net of the MPMS magnetometer during the field ramp. 
The supercurrent density evaluated from Eql2] (shown by 
stars in Fig. [8|) compares well with the magnetization 
data. Finally, the inset in Fig. [5] compares the supercur- 
rent obtained from magnetic measurements with direct 
transport measurements. There is good agreement be- 
tween all datasets. We therefore conclude that, similar 
to NdFeAsOi-ajFj^^i, the supercurrent density in single 
crystals Ba(Feo.93Coo.o7)2As2 is relatively low even at 
low temperatures, j (T = 5K, = 0) ~ 2.6 ± 0.2 x 10^ 
A/cm^, which is comparable to Y-Ba-Cu-0 single crys- 
tals. As for the functional form of jc{T), we can use 
the generalized power law that is usually used for su- 
perconductors, j(T) = j{0) (1 — (T/Tc)^)", which gives 
the best fit to the experimental data (most data points 
- solid symbols) with j(0) w 0.88 MA/cni^, p = 1/4 and 
p = 1. An attempt to fit with p = 1 and n = 3/2 is also 
shown in Fig. [S] and is inferior to the former fit, which 
is reflected in an order of magnitude larger adjusted 
value. Next we present unexpected results obtained in 
a magnetic field. Using the annealing method described 
above for the remanent state, but this time annealing the 
induced critical state in finite fields, we have obtained a 
set of j [H — const, T) curves for different values of H . 
The result is shown in Fig. O There is an apparent pro- 
nounced decrease in the current density marked by open 
circles. Similarly, the temperature at which jab — » is 
also marked. While the meaning of the former signature 
will become clear below when we discuss the — T di- 
agram, the latter feature is obviously the experimental 
irreversibility line above which pinning is negligible and 
the vortex liquid sets in. 
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FIG. 9: (Color online) j (T) measured in a single crystal 
Ba(Feo.93Coo.o7)2As2 upon warming up at different values of 
the magnetic applied field as explained in the text. Open cir- 
cles mark the temperature of the vanishing current as well as 
the temperatures at which j (T) curves exhibit an apparent 
change of behavior. 



FIG. 10: (Color online)Time - dependent magnetization mea- 
sured for 22 hours at T = 15 K at ii" = 1 T after ramping 
field from 5 T. Lower inset: fit to the collective creep model, 
Eq. |4l with parameters indicated. Upper inset: In — In scale 
graph of relaxation curves obtained at 1 T and 3 T. 



B. Magnetic relaxation 

When the fishtail effect is observed, one must examine 
the magnetic relaxation of the critical state as a function 
of temperature and magnetic field. While some theories 
suggest a static explanation for the fishtail effect, which 
most likely works in alloys, ceramics, and materials with 
secondary phases, in single crystals it seems that a dy- 
namic scenario is more plausible, at least in the high-Tc 
cuprates.— 1^ Within collective pinning and creep models, 
the magnetic relaxation rate decreases with the increase 
of a magnetic field. Since the true critical current den- 
sity, jc{H), also decreases with the increase of a mag- 
netic field, an apparent fishtail develops in the measured 
magnetizationiS Therefore, the fishtail phenomenon is 
a direct consequence of the collective pining and creep 
theory. However, more detailed measurements above the 
peak position have forced us to revisit this picture and 
use a refined model that involves both field-dependent 
relaxation rates and a transition from a collective flux 
creep mechanism at the lower fields to a plastic creep 
(mediated by dislocations in the vortex lattice).'*^ The 
latter predicts an increasing relaxation rate with the in- 
crease of the magnetic field, whereas the former predicts 
a monotonic decrease of the relaxation rate. 

Based on the various models of magnetic relaxation, 
Griessen et al. has suggested a very useful generic form 
of the barrier for flux creep.~ 



U{3) 



Uo 
a 



1 



(3) 



tional forms of U {j / jc) if the exponent a is allowed to 
have both negative and positive values. For a — ~1, Eq. 
[3] describes the classical Anderson-Kim barrier, which is 
simply linear in j/j<r^-^- For a = —1/2 the barrier 
for plastic creep^ is obtained. Positive a describes col- 
lective creep barriers^. In the limit a ^ 0, this for- 
mula reproduces exactly the logarithmic barrier.*^ An 
activation energy, written in the form of Eq. [3l results 
in a so-called interpolation formula for flux creep^ ob- 
tained when the logarithmic solution of the creep equa- 
tion U (j/jc) = Tln(l + t/tof^ is applied (for a ^ 0), 



M(0 = M.(l + |^ln(l 



(4) 



This formula describes all the widely-accepted func- 



For a = 0, a power-law decay is obtained: M (t) = 
Mc{to/t)" , where n — T/Uq. A generalization of the 
flux creep theory beyond the logarithmic solution with 
applications to collective creep and the fishtail effect are 
given by Burlachkov, Giller, and Prozorov.-ii 

Indeed, no single equation for relaxation is applica- 
ble for all values of j/jc- Within the collective creep 
theory, the magnetic relaxation has different functional 
forms (roughly described by using different Uq and a in 
EqlSj depending on whether one considers the regime of a 
single- vortex, small vortex bundles, or large ones^. Try- 
ing to find the best fit to Eq.|4]in the widest time- interval, 
we have determined that the small-bundle regime de- 
scribes the data shown in Fig. [TO] quite well. The fit 
is shown in the lower inset. The fit parameters are rea- 
sonable within the collective creep model, as a = 5/2 
(fixed during fitting), Uq = 564 K, and to = 2 fis. 

We now examine relaxation magnetic fields above and 
below the peak position, Hp{T). As can be seen in the 
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FIG. 11: (Color online) Magnetic relaxation measured in a 
single crystal Ba(Feo.93Coo.o7)2As2 at different values of the 
applied magnetic field shown along with the M{H) loop mea- 
sured at the same temperature, T = 15 K. Right axis shows 
the logarithmic relaxation rate as function of an applied field. 
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FIG. 12: (Color online) Comparison of the relaxation rate, R, 
as function of an applied magnetic field measured in a single 
crystal Ba(Feo.93Coo.o7)2 As2 in two different orientations at 
T = 15 K. 



upper inset of Fig. [TOl already just above the peak, the 
relaxation rate is faster than before the peak. Its func- 
tional form is also quite different. The magnetic field de- 
pendence on the relaxation rate is presented in Fig. [TT] 
Here, each M [H = const, t) trace is plotted along with 
the regular magnetization loop. 

A useful characteristic of the flux creep is the logarith- 
mic relaxation rate, R = — \d\nM/d\nt\. From Eq[4]we 
obtain. 



trend when a magnetic field is applied along the ab— 
plane. We therefore conclude that the collective creep 
model fails above Hp. 

Another difference between the relaxation rates comes 
from examining the temperature dependence of R. Fig- 
urefTSl compares R (T) measured in zero field (self-field of 
the remanent state) for three different orientations (along 
three of the edges of the crystal) . When vortices are par- 
allel to the ab— plane, the relaxation rates are practically 



(5) 



Even at the longest time that we have probed (22 
hours), the logarithmic term is equal to 1.6, so the relax- 
ation rate only decreases by a factor of about 2.6 com- 
pared to the initial (but not too close to t = 0) rate, 
R ~ T/Uo, which at 15 K and 1 T is about R = 0.027. 
This value is much larger than the rate for conventional 
superconductors, but is quite similar to high- Tc com- 
poundsj^ Thus, single crystals of Ba(Feo. 93000.07)2^82 
exhibit "giant flux creep". The magnetic field depen- 
dence of R is summarized in Fig. [Tl] (right axis). At 
first, it decreases with the increase of H, up to the field. 
Hp, corresponding to the maximum in the fishtail M (H) 
curve. Above Hp the relaxation rate increases with the 
increase of H, which contradicts the collective creep the- 
ory. To clarify whether this behavior is related to the 
fishtail phenomenon, we have measured the magnetic re- 
laxation rate with the magnetic field applied along the fl- 
edge and we did not observe a fishtail feature (Fig. [3]). 

Figure [T^ compares the relaxation rates in the two ori- 
entations. While above Hp, R{H) increases with H for 
measurements along the c— axis, it shows the opposite 
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FIG. 13: (Color online) Logarithmic relaxation rate as 
function of temperature measured in a single crystal 
Ba(Feo.93Coo.o7)2As2 in the remanent state, ii" = 0, in three 
orientations of the trapped vortex direction with respect to 
the sample. For comparison, R (T) at = 3 T is shown. 
(Note different origins and scales of the left and right axes). 
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the same. For magnetic relaxation in the remanent state 
obtained along the c— axis, R{T) shows a similar func- 
tional dependence and is shifted slightly to lower temper- 
atures. In contrast, the relaxation rate obtained above 
Hp shows a very different behavior - increasing with the 
increase of temperature. The bell-shaped R (T) are ob- 
served in many superconductors with weak pinning. No 
straightforward explanation exists, but it seems plausible 
that while the bare barrier for the magnetic relaxation, 
UT.j=o, decreases with the increase of temperature, the 
ratio j / jc at which the relaxation takes place in the ex- 
periment with the fixed time window Ai decreases with 
the increase of temperature. Therefore, the actual bar- 
rier experienced by the vortex bundles, Eq. O may be- 
come non-monotonic with temperature. In other words, 
the exponent a is temperature dependent and from the 
collective creep theory it follows that it is nonmono- 
tonic with the bundle size, peaking at a — b/2 for small 
bundles^. By substituting this into Eq.Owith the appro- 
priate UT.j=Q (T), the nonmonotonic R (T) can be repro- 
duced. This provides additional indirect evidence for the 
applicability of the collective creep approach to vortex 
relaxation in single crystalline Ba(Feo.93Coo.o7)2As2 at 
fields below Hp. Using the same arguments, one can show 
that plastic creep, which corresponds to a = —1/2 can- 
not lead to such behavior and the relaxation rate should 
just increase with temperature, which is what we observe 
in Fig. M 

The idea of plastic vortex creep was introduced when 
a similar failure of the collective creep model was found 
in YBCO^ and the electron-doped cuprate supercon- 
ductor, Ndo.85Ceo.i5Cu04_2; (NCCO)^, single crystals. 
Similar measurements were used to extract E (j) char- 
acteristics without attaching contacts and confirmed the 
plastic creep model^. These and later works showed 
that in superconducting crystals with relatively weak pin- 
ning, the crossover from elastic to plastic creep always 
accompanies fishtail magnetization. However, it should 
be noted that the fishtail feature itself is still due to pecu- 
liarities of the collective creep when the relaxation rate is 
faster for lower fields, H < Hp. Within the plastic creep 
scenario, thermally activated vortex motion is not due 
to jumps of vortex bundles, but due to sliding of disloca- 
tions in the vortex lattice (even in a very disordered state, 
one can consider a local ordered arrangement of vortices 
and dislocations as primary defects). The main result 
of this approach is that the barrier for magnetic relax- 
ation in the plastic creep regime remains finite as j — > 0, 

Upi (BJ/jc) = Upi (B) (l - ^JThc^ and decreases with 

increasing magnetic field as Upi {B) ^ 1/^/B, because 
the elementary movement of a dislocation is proportional 
to the intervortex distance, a ~ 1/^/B. Collective creep 
barriers, on the other hand, increase with increasing mag- 
netic field as Uc {B) ^ B'^ , where ly is positive and de- 
pends on the particular pinning regime. The j depen- 
dence is given by Eq|3]with a positive a that results in a 
diverging barrier, Uc {] I jc) ~ {j I jc) ■ At low fields, the 




T(K) 

FIG. 14: (Color online) Vortex H — T phase diagram obtained 
from various features in the described measurements and the 
solid lines are the fits as discussed in the text. 



collective barrier is always smaller than the plastic barrier 
and this channel dominates the vortex dynamics. How- 
ever, as the magnetic field is increased, a crossover from a 
collective to a plastic channel of vortex relaxation occurs 
when Upi {B,j/jc) becomes smaller than Ucoi {B,j/jc)- 

Finally, Fig. [T3] presents the vortex phase diagram of 
Ba(Feo.g3Coo.o7)2As2 single crystals compiled from our 
results. It closely resembles diagrams for the cuprate su- 
perconductors, in particular YBCO single crystals. The 
prominent features are the irreversibility line that is quite 
distant from the Hc2 (T) line and a clear crossover line 
between collective and plastic creep regimes. As is usu- 
ally the case, the data can be described by the general 
power law, H,,{T) = i?,,(0) (1 - (T/TJ^)", shown in 
Fig. [13] by solid lines. We obtained for the fishtail peak 
position, Hp{T), Hp{0) « 13.6 T, n = 3/2, and p = 1. 
For the irreversibihty fine, Hirr{T), Hirr{0) ~ 15.2 T, 
n — 3/2, and p = 2. For the second critical field, Hc2{T), 
-ffc2(0) ~ 66.7 T, 71 4/3, and p = 1. 



IV. CONCLUSIONS 

In conclusion, detailed measurements of global and 
local electromagnetic properties clearly show that the 
vortex behavior in single crystalline superconducting 
Ba(Feo.93Coo.o7)2As2 is similar to that found in the high 
Tc cuprates. In particular, we report the presence of 
a fishtail feature in AI [H, T ~ const) loops and find its 
signature in M [H — const, T) measurements. Further- 
more, magnetic relaxation measurements are consistent 
with the collective pinning and creep models (weak pin- 
ning and fast creep) and suggest a crossover into the plas- 
tic creep regime in fields exceeding the value correspond- 
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